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Abstract 

After the discovery of the phenomena of hght-induced excited spin 
state trapping (LIESST), the functional properties of metal complexes 
have been studied intensively. Among them, cooperative phenomena 
involving low spin-high spin (spin-crossover) transition and magnetic 
ordering have attracted interests, and it has become necessary to for- 
mulate a unified description of both phenomena. In this work, we 
propose a model in which they can be treated simultaneously by ex- 
tending the Wajnflasz-Pick model including a magnetic interaction. 
We found that this new model is equivalent to Blume-Emery-Griffiths 
(BEG) Hamiltonian with degenerate levels. This model provides a uni- 
fied description of the thermodynamic properties associated with var- 
ious types of systems, such as spin-crossover (SC) solids and Prussian 
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blue analogues (PBA). Here, the high spin fraction and the magneti- 
zation are the order parameters describing the cooperative phenomena 
of the model. We present several typical temperature dependences of 
the order parameters and we determine the phase diagram of the sys- 
tem using the mean-field theory and Monte Carlo simulations. We 
found that the magnetic interaction drives the SC transition leading 
to re-entrant magnetic and first-order SC transitions. 
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I. INTRODUCTION 



Recently the functionality of metal complexes attracts much interest. In par- 
ticular, the cooperative phenomena where both spin-crossover (SC) and magnetic 
ordering are involved have been studied intensively. Molecular switching of inor- 
ganic solids is a typical issue of the vibronic lability of molecular units, originally 
introduced in chemistry by J. A. Ammeter [1,2]. 

The concepts of molecular bistability were given by O. Kahn [2], and recent 
reviews on this subject can be found in the references [3-6]. In addition, due to the 
characteristic of reversible control of the magnetic properties by light, temperature, 
pressure, and magnetic field, SC and PBA compounds are considered to be promising 
candidates for information storage. [7-9]. 

system low-temperature sate High temperature State ref. 

SC system Fe" {3d^) Fe^s (S=0) [t%] Fe'^s (S=2) [4^e%] [10-12] 

PBA Fe,Co Fe[^s (S=0) Co[^^(S=0) Fe}^|(S=l/2) Co}^s(S=3/2) [13-15] 

In SC, the electronic structure changes between the high spin state (HS) and 
and the low spin state (LS) in an atom which causes changes of magnetic property 
and also volume of the molecule. On the other hand, in PBA the charge transfer be- 
tween Fe and Co causes a LS-HS (spin-crossover) transition of Co, and the magnetic 
moments of the Co and Fe couple antiferromagnetically in the high spin state (HS), 
which results in local ferrimagnetic moment. Although the LS-HS transition of PBA 
is different from simple SC transitions, we can regard PBA as a SC system in the 
broad sense. A new aspect of PBA is the magnetic ordering at low temperatures, 
where the magnetic interaction between spins is important. 

As another type of phase transition due to the charge transfer, other 
new inorganic systems, showing magnetic changes after the spin transi- 
tion, have been reported in literatures, such as mixed-valence iron complex 
(nC3H7)4N[Fe"Fe"^(dto)3](dto=C202S2) [16,17], where the charge transfer causes 
a change of degeneracy, and a phase transition occurs between a high temperature 
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and low temperature states. In this system, a magnetic transition is also observed at 
a low temperature. There, the spin-crossover of each atom does not occur. However, 
the entropy effect plays an important role in the structure change associated with 
the electronic configuration. It is an interesting problem to study the mechanism of 
combined phenomena of the structure change and magnetic transition. 

In SC transition, a cooperative distortion of the lattice occurs, which is associated 
with a large change in the average metal-ligand distance. Indeed it increases about 
10% at the molecular level in the conversion from the LS to HS state which causes 
a volume expansion of the molecules. Thus, the SC phenomenon at the microscopic 
scale is due to the coupling between the electronic and the vibrational structures 
[18] of the molecules. That is the intra-molecular vibronic couphng. The couphng 
between the lattice distortion and the change of the electronic states of the molecules 
induces the cooperativity of SC sohds. This is the microscopic origin of the elastic 
interaction which is responsible for the SC transition. Prom the experimental point 
of view, optical and magnetic properties change during the SC transition, so that the 
switching properties can be followed by magnetic or optical (absorption, reflectivity) 
techniques. 

For the analysis of these structural transition, it is largely admitted that the 
bistable properties at the molecular level are adequately described through a molec- 
ular configurational diagram, i.e. a plot of adiabatic energies versus distortion co- 
ordinate of the molecular system. In Fig.l, we show a unified three-states config- 
urational diagram suited to SC and PBA in the case of the low spin (LS) ground 
state. It is noteworthy that the effect of environment in molecular sohds affects 
the configurational diagram. For example, an external pressure mainly increases 
the £'(HS) — E{LS) energy gap, and thus reduces energy barrier between HS and 
LS states in SC solids. Consequently, it raises the transition temperature (Teq) and 
it decreases the lifetime of the metastable HS state [19-21]. In the case of PBA, 
the same effect of pressure is also expected on the spin transition of cobalt [22], 
which would cause variation of the magnetic properties under high pressures. In the 
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mixed- valence iron complex (nC3H7)4N[Fe"Fe"^(dto)3](dto=C202S2), the pressure 
effect has been reported [23]. 
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FIGURES 



FIG. 1. A unified three-states configurational diagram. The non-degenerate ground 

state is of low spin and the excited degenerate states are of higli spin. 

To describe SC systems, a two-states Ising-like model proposed by Wajnflasz and 
Pick (WP) [24] is currently used. In this phenomenological model, the HS and LS 
microscopic states of the molecule are associated, respectively, with the eigenvalues 
-|-1 and —1 of the fictitious spin operator s with their respective degeneracies g+ and 
g^. Here, we denote hy g — ^ > 1 the degeneracy ratio, by A > the ligand field 
energy (such that the LS state is the ground state at OK), and by J' the interaction 
between the molecules. J' is assumed to have an independent value on the spin 
states. Taking into account the isomorphism [25] between the Ising model with 
degenerate levels and an Ising model under a temperature dependent "field" [26], 
the WP Hamiltonian is written as 



where J2 is the sum over the interacting neighbors and s takes 1 (HS) or -1 (LS). 



J' > is the interaction parameter between the SC units. 

As we mentioned before, the interaction in SC solids originates from elastic 
strength [27]. However, in the WP model it is expressed through the phenomeno- 
logical term —J'J2 SiSj. Here J' is not an exchange magnetic coupling. This term 



expresses the situation that an attractive force (i.e. a ferroelastic interaction) acts 
between neighboring sites in the same state (LS-LS or HS-HS) and a repulsive force 
acts between those in different states (LS-HS), which leads to bistability of the LS 
and HS states. A more realistic treatment which includes a spin-phonon interaction 
leads to obtain a similar Hamiltonian as (1) with effective parameters J' and A 
which depend on the distortions of the molecules at sites i and j [28] . Thus, it is 
important to consider here the parameters J' and A as effective parameters which 
have elastic origin, even if they are taken as constants in the model in order to make 
the present treatment easy. 




(1) 
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The degeneracy ratio (g) between the HS and LS states is related to the molar 
entropy change (AS') due to the total spin-conversion which is given by AS — 
NkB In g , where /cb is the Boltzmann constant and N the Avogadro's number. 
Because of both electronic, vibrational (intramolecular) and phonon factors, the 
experimental entropy changes by AS = 30 — 70JK~^mol~^ [4]. This corresponds 
to the following degeneracy values g = 36 — 4500, which indicates that g ^ 1. 
Generally, g depends on temperature through the temperature dependence of the 
phonon density of the lattice [29]. Here, we simply take g as a constant. 

The ratio of molecules in the high spin state, called HS fraction, is denoted by 
riHs, which is expressed as a function of the "fictitious magnetization" < s > such 
as TT-Hs — (< s > +1) /2. Despite of the drastic simplification of the realistic vibronic 
levels scheme by the two-level system, this model permits to describe with amazing 
success most of the quasi-static properties of the SC solids. For example, it explained 
consistently the change from a smooth transition to a first order transition [24,26,30] 
of spin conversion between the LS and HS states as a function of the parameters of 
the model. It also well explained the two-step spin crossover behavior [31-33]. 

In order to study the magnetic properties in PEA sohds, etc., we need to go 
beyond the 2-states model. In previous studies [34,35], focusing on the condition of 
the reversible switching between nonmagnetic LS state and the magnetic metastable 
ordered HS state which was observed in PBA [13], the following three state model 
called Blume-Capel model [36,37] was investigated: 

H = -JJ2S^S, + AJ2Sl (2) 

In this model, the spin operator Si has three states: —1, 0, +1, where (±1) is asso- 
ciated with the LS nonmagnetic (HS magnetic) state. Here, J is really an exchange 
magnetic coupling and A plays a role of a ligand field. Reversible switching was 
demonstrated by a suitable change of parameters of the model assuming photon's 
effect is renormalized into the parameters. In addition, from a dynamical viewpoint, 
it has been pointed out that existence of multi-time scale dynamics is essential for 
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the switching in both directions (magnetic <-> diamagnetic) . That is, the difference 
between the time scales of the spin flip process and of structural change (LS'S-^HS) 
process is important. This concept is very important to investigate photomagnetism 
of PBA [38]. 

To reproduce the temperature-induced SC properties of PBA, the information of 
the difference in degeneracies between the high temperature and the low temperature 
phases and the nonmagnetic interaction between molecules (which is originated from 
the lattice) have to be taken into account. This is the motivation of this paper. 
Here, we propose a model taking into account both mechanisms of the temperature- 
induced SC transition and of the magnetic transition. Then, we extend the 3- 
states model [34,35] (given by the Hamiltonian (2)) to include the mechanism in the 
Wajnflasz and Pick model (Hamiltonian (1)). As we will see in the next section, 
this present model can be regarded as an extended Blume-Emery-Griffiths (BEG) 
Hamiltonian [39]. 

It is worth noting that for the elastic interaction between molecules, alternative 
scenarios to describe the SC phenomenon have been proposed recently. Some of them 
are based on pseudo-spins coupled by springs [28] . The originality of this approach is 
that the elastic constant between two neighboring atoms depends on their electronic 
states. Indeed, due to the expansion of the lattice in the HS state, the elastic force 
constants are weaker in the HS state than in the LS state. This approach leads to 
a Hamiltonian which is equivalent to an Ising-like model in which the interaction 
parameter depends on the local distortions of the molecules. An extension of this 
scenario including the intramolecular potential energy of the molecules has been 
published recently by one of the authors [40]. However, due to the complicated 
structure of the Hamiltonians, these models have been investigated only in the one 
dimensional case, but we expect new development of this kind of descriptions in 
near future. 

This paper is organized as follows: in section II we introduce the model, section 
HI is devoted to study the equilibrium properties obtained by a mean- field theory. In 
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section IV, results of Monte Carlo simulations are presented. In section V, we discuss 
the obtained results in comparison with available experiments in the literature and 
we conclude. 



II. THE HAMILTONIAN 

Here we adopt a model based on a three states Ising-hke Hamiltonian oi S — 
0, +1,— 1. The state S = denotes the non-magnetic state, which corresponds to 
the electronic configuration FcLgCoLg in the case of PBA, and simply Fcls f*^^ SC 
systems. The states S = ±1 are associated with magnetic states. The magnetic state 
corresponds to the ferrimagnetic state FeLgCo^g for PBA compounds, and Fe^g in 
the case of SC sohds. Due to large expansion of the coordination sphere of the 
molecule during the spin transition, the vibrational structure [4] and the electronic 
properties between the HS and LS states are strongly affected. Consequently, the 
degeneracies of the non- magnetic state and the magnetic state are different. Let us 
denote by u and r the degeneracies of the states S — and 5" = ±1, respectively, and 
by g the ratio ^. Here we assume that the degeneracies of the states with S = ±1 
are the same and g is temperature independent. 

Because the high and low spin states are described by -S" = ±1 and -S" = 0, 
respectively. The Wajnflasz and Pick model (1) in the present three states approach, 
in which the degeneracies are included as entropic terms [24,26,25,41], is given by 

i^wP - -J'T^i^Sl - 1)(25^ - 1) + (a - ^In^) E(2^^ - !)■ (3) 

Adding the magnetic interaction —JJ2i^jSiSj and the Zeemann energy —hJ2iSi, 
we obtain the general Hamiltonian: 

H = - J^^SiSj - hY^Si - KY, S^S] + {D- k^T In g)Y^S^ + const. (4) 

where K = 4 J' (>0) and D = 2 A + 4zJ' >0 {z : the coordination number). Here 
only the nearest neighbor interactions are considered for the pairing i ^ j. The 
magnetic interaction between the spins is only relevant when both spins Si and Sj 
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are in the magnetic state, i.e. 5" = ±1. If we do not consider the degeneracies, i.e., 
u — r — 1 and g = 1, this Hamiltonian corresponds to the well known BEG model. 
[39] Here, we adopt the form of Eq. (4) to refer to past works on BEG model. 
Indeed, BEG Hamiltonian was extensively studied to describe the thermodynamical 
properties of the tri-critical point in He^-He^ mixtures [39]. Later, Sivardiere and 
Lazerowitch [42] used an extended version to lattice gas model of BEG Hamiltonian 
to study the condensation and phase separation in binary fluids. It is also interesting 
to mention that even recently, BEG Hamiltonian was also apphed [43] to describe the 
ferroelectric order thermally induced in donor-acceptor materials like TTF-TCNQ 
[44]. This kind of materials are also intensively studied for their photo-induced 
properties [45]. 



III. EQUILIBRIUM PROPERTIES 

In experiments, so far the spin crossover and magnetic ordering in PBA take 
place in the range of temperature 10 —300 K. As typical values of the parameter 
of the model, extracted from the experiments [15,46], we can refer to the exchange 
interaction zJ ~ 20K , the hgand-field energy gap D 500 K, the elastic interaction 
zK is about 100— 200K and the degeneracy ratio g ~ 100—200. Here z is the number 
of the nearest neighbor pairs. The present model describes the systems of SC and 
PBA solids as special cases: 

(i) for SC systems, the magnetic interaction between the magnetic HS units is 
negligible because of the large distance between them (pa 10 A) [47]. 

(ii) In the case of PBA, all the terms of the model should be taken into account, 
because they are effectively present. However, as far as we concern the equilibrium 
thermodynamics, only the electron-transfer contributes to the thermal bistability 
between the states FclsCols and FclsCohs- Indeed around the spin crossover tran- 
sition temperature, experimentally around 200— 300K, there is no magnetic ordering. 
This fact suggests that the magnetic coupling zJ, which is about 20K, is very small 
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in comparison to the ligand-field energy gap D 500 K. However, in the near fu- 
ture, because of an amazing developments of synthesis of materials in this subject, 
we may expect new systems in which the magnetic and SC properties could com- 
pete. Thus, we study a general structure of the phase diagram caused by K and J 
in the present paper. 



A. The SC case 

As we mentioned above, in the case of SC, the magnetic interaction is irrelevant 
and the Hamiltonian (4) without the magnetic field and magnetic interaction can 
be written in the following form 

Heg = -K^^S^S] + {D- k^T\ng)Y,Sl (5) 

It is straightforward to see that a flipping of any spin {Si —>■ —Si) does not change 
the Hamiltonian. Therefore, as expected, the magnetization m —< S > should be 
zero at any temperature. Introducing a new variable Uj = 2Sf — 1, the Hamiltonian 
5 reduces to WP model [24], which leads to a first order phase transition as it is 
well known from the literature [26,30,48,49]: 

Ha^ -^2^ (Ti(Tj H ^ cr^ + const. (6) 

i^i i 

In a mean-field approximation, the order parameter of the system is q —< S'^ >. 
It corresponds exactly to the HS fraction of molecules in the system. The self- 
consistent equation accounting for the thermally driven first-order transition be- 
tween the macroscopic LS {q ~ 0) and HS {q ~ 1) states, is given by 

D-zKq 

So-called transition temperature Teq is obtained by setting q — 1/2 in the equa- 
tion (7), or more simply by setting the effective field D — kBTln2g — zK/2 of the 
Hamiltonian (6) equal to zero. Thus we have 
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, _ D~zKl2 

It is interesting to notice that occurrence of the first-order transition is only 
possible when the condition Tcq < Tj, where Ti denotes the critical temperature of 
the pure spin 1/2 Ising model in the same lattice. In the mean-field approximation, 
the Curie temperature is given by Ti — zK/Ak-Q. 

Two types of behavior of the order parameter can be obtained in the mean-field 
approach depending on the ligand-field energy gap D and the interaction K: 

(i) a first-order spin-crossover transition, with a hysteresis loop, when zK > 
4D/{2 + \n2gy, 

(ii) a smooth conversion between the two states when zK < AD/ {2 -\- ln2gf). 
The temperature dependence of q is summarized in Fig. 2. 

FIG. 2. Temperature dependence of the high spin fraction q for values of the energy 
gap D. From the left to the right, D = 90K, lOOK, 120K and 150K with the fixed values 
of parameters K = 90, and g = 10. Smooth (first-order) transitions are obtained for 
D > lOOK {D < lOOK). The part of the curves with negative slope dq/dt represents 
unstable states. 



B. General case. 

Here, we consider the general model of SC solids in which the magnetic moments 
of the HS state interact magnetically, which may cause a magnetic order. Therefore, 
in this section we consider the Hamiltonian (4) with various parameters. The model 
imposes two order parameters m —< S > and q —< S'^ >, respectively, associated 
with the magnetization and the high spin fraction in the system, li K ^ J as the 
case of PBA, we can deal with a multi-scale energy problem in the present frame 
work: high energy scale for the spin-crossover phenomena (at high temperature), 
low energy scale for magnetic phenomena (at low temperature). 

Now, let us study properties of the system (4) by a mean-filed theory (MFT) 
which gives the phase diagram of the model. The one site Hamiltonian for this 
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system is given by 

Hi^- (zJm + h)Si- {zKq - D + UbT In g) Sf. (9) 

The associated mean-field free energy per site is given by 
F{m, q, T) = -zJm^ + -zKq^ - kj^T In (l + g/^C^-^^-D+teT in 25) ^^^^ ^ ^ _ 

(10) 

Minimizing the free energy (10) with respect to the order parameters m and q, we 
easily obtain the following coupled self-consistent equations on m and g, 

= — T^n? — n_Lt. o ^ r^TTi ;^ — - = qtsjahp{h, + zJm) (11) 

and 

cosh /3 (/i + zJm) 

^ = e-/3(zifq-I3+feBTln2s) + cOsh /5 (/l + zJui) ' ^ ' 

li J — h — 0, we reproduce naturally the case of SC previously studied. Equa- 
tions (9,11,12) have the same structure as that found in BEG model. However, due 
to the temperature dependent ligand-field (anisotropy in BEG model) we expect 
new behavior for the order parameters m and q and also a different phase diagram. 
Putting g — 1 gives exactly the same self-consistent equations as that published by 
Mukamel and Blume [50] in their study of the BEG model to describe the thermo- 
dynamical properties of the tricritical point in He^-He^ mixtures. Our aim here is 
to analyze the effect of the magnetic interaction on the spin-crossover phenomena 
{g ^ 1). Therefore, wc study the following two cases: 

1. J ^0,h^O,K y^O,D^O; 

2. J y^O,h^O,K y^O,D y^O; 



1. J = 0,hj^O,K j^O,D j^O case 



This case corresponds to the situation of the SC system under a magnetic field 
[51,52]; the spins are too far to interact each other. In such situation, the system is 
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paramagnetic or diamagnetic, and the magnetization has no strong effect. However, 
we study the effect of the magnetic field on the transition temperature and on 
the hysteresis loop in detail. In the present situation, the self-consistent equations 
(11,12) become 

m — qtanhph (13) 

and 

cosh/5/i , , 

^ = e-p{zKq-D+k^Tln2g) + cOsh/3/l' ^ ^ 

It is clear from Eqns. (13 and 14) that the magnetization depends linearly on the 
HS fraction in the system. In the limiting case, where the magnetic field is very small 
in comparison with the transition temperature, the magnetization is m ~ q(3h. In 
the following, we will demonstrate that the magnetic field stabihzes HS (magnetic) 
state and decreases the transition temperature T^q of the system. So, we expand 
the equation (14) around the transition temperature of the first-order transition, for 
which q — 1/2, by considering that ph <^ 1. After some simple calculations, we 
arrive at the following expression of the field dependence of transition temperature 

Teq(/l) 

Te,(/.) = Te,(0) - (15) 

Here, Teq(O) = ^f^l^2g transition temperature of the system under zero applied 

magnetic field. The shift of the transition temperature 

AT = - = - fl6) 

'^'^'^ feTeq(0)ln2^? AH{0) 

depends on the ratio between the square of the magnetic field and the enthalpy 
change, Aif(O), at the transition. These results are in good agreement with available 
experimental results and the microscopic analysis on SC [52,53] phenomena. It 
follows from this result that, lowering the transition temperature Teq{0) enhances 
the magnetic field effect on SC systems. 
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FIG. 3. Effect of the 'magnetic' field h {h = OK, h = 40K) for K = 200K, D = 300K, 
g = 150) on the temperature dependence of the HS fraction q{T) in the case of the 
first order phase transition. The magnetic field stabilizes the HS phase and polarizes the 
magnetic system in which it induces a nonzero magnetization (dashed line). 

On a hysteresis loop, see Fig. 3, the magnetic field has as effect to shift the 
whole hysteresis to the low temperature side, by acting differently on the ascending 
and descending metastable branches. Denoting, by (TJ^) and {Tq), the 
transition temperatures under the field h {h — 0), at the ascending and descending 
branches, respectively, one can deduce approximately from the equation (15) that 
the associated shifts are given 

feAT+ = -tJ^i^ and /cbAT" = -jfj^- (17) 

Since Tq^ > Tq , it follows that AT^ > AT"*", which indicates that the magnetic field 
causes a non-symmetric distortion of the hysteresis loop of SC solids by stabilizing 
the HS and destabilizing the LS state. 

To give an order of magnitude of the shift of the hysteresis loop, we consider 
the case of Co-'^-'^(H2(fsa)2en)(py)2 [H2(fsa)en = 

N,N'-ethylene bis(3-carboxysalicylaldimine,py) = pyridine] spin-crossover solid [54] 
which shows an abrupt thermally induced spin- crossover transition centered at 126K. 
For this system, the magnetic spin moments corresponding to the HS and the LS 
states are 3^ — 3/2 and Si^ — 1/2, and the energy gap is 700 K [55]. In our Hamil- 
tonian, the parameter h denotes h = gi^fiBS, where /x is the Bohr magneton, and B 
the applied magnetic field. Here gi, is the Lande factor g'L = 2 for both LS and HS 
states. For Co^^(H2(fsa)2en)(py)2, we found that the shift ATeq ^ —4.9 K under a 
static magnetic field of 32 Tesla. 
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2. J ^0,h = 0,K ^0,D j^O case. 



In this case, the HS units created thermally or optically interact through the ex- 
change coupling J. We study effects of this magnetic interaction on the charge trans- 
fer involving SC transition, under zero applied magnetic field. In MFT, the thermo- 
dynamical properties are obtained by solving the following coupled self-consistent 
equations: 

m — qtanh (3zJm (18) 

and 

cosh [3zJm 

^ = g-/3(^i^q-i?+feBTln23) + cosh /3^Jm' ^^^^ 

Here, it is clear from the equation (18) that the system may order magnetically, 
depending on the thermal behavior of the order parameter q. Indeed, at Curie 
temperature Tm, we have 

q{TM)PMzJ = 1. (20) 

It is interesting to notice that the thermal behavior of this complicated situation 
may be discussed through the comparison of the transition temperatures of the 
magnetic and spin- crossover systems. At Curie temperature Tm, the magnetization 
777.(7) shows a second-order phase transition, and the HS fraction q{T) exhibits a 
singularity. 

Below Tm spontaneous magnetization appears where the system is globally in the 
HS phase. As the temperature decreases, the HS fraction decreases smoothly and the 
magnetization disappears at Tc(< Tm) with a second-order transition for small value 
of J where q{T) exhibits a singularity. When J is large, the magnetization remains 
until Teq and there disappears discontinuously at the first order phase transition 
temperature. 
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FIG. 4. Temperature dependence of m{T) (thin curves) and q{T) (bold and dashed 
curves) in the case zK = 90K, D = 150K and g = 10 for the following values of the 
magnetic interaction. From the bottom to the top: zJ = 48K, zJ = 57.7K, zJ = 58K, 
zJ = 59K, zJ = 61K, and zJ = 63K. The solid bold curve denotes the HS fraction q for 
J = OK, in which m = 0. Note that the magnetization m{T) is equal to zero in the case 
zJ = ASK. 

When the SC and magnetic transition temperatures are of the same order of 
magnitude, the magnetic system (i.e., m{T)) perturbs the "lattice" properties (i.e., 
q(T)) as seen in Fig. 4, where the bold solid hne denotes the temperature dependence 
of the order parameter q for J = 0. There all the simulations have been performed 
varying the magnetic coupling J at fixed values of K and D: zK = 90 K and 
D — 150K. For small values of J, as zJ — 48K, no magnetization appears and only 
smooth SC transition occurs. When zJ has values between 57.7 and 59 for example, 
we observe a second-order re-entrant phase transition on the magnetization m with 
the transition temperatures Tc and Tm- The associated HS fraction shows singularity 
at each transition temperature. Increasing zJ, such as zJ = 63K, the magnetization 
presents a first-order transition at Tc and a second-order transition at Tm. In such 
case, the HS fraction exhibits a first-order transition at Tc and a singularity at Tm. 
Thus, it appears clearly that the magnetic coupling causes the first-order transition 
in the elastic properties (i.e., in q{T)). In the following we discuss the three regions 
identified in Fig. 4: 

(1) a high temperature region in which g ~ 1, m = 0, corresponding to a 
paramagnetic HS state; 

(2) an intermediate region Tq < T < Tm, where > 1/2, m 7^ 0, indicating the 
existence of a predominant HS state with a ferromagnetic order; 

(3) at low temperature T < Tc, we have a phase of LS. There, we find two sub- 
regions. At relative high temperature region, q is still large, where many high spins 
remain, which are made of isolated magnetic states 5" = ±1 in a sea of diamagnetic 
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states S — where the system is paramagnetic. On the other hand, at low temper- 
ature q is very small and the system is almost the perfect diamagnetic state. In this 
low temperature phase in which m = and q < 1/2, we obtain a reentrant magnetic 
phase. Between the regions of (2) and (3), we find a first order SC phase transition 
driven by the magnetic interaction. Due to the magnetic interaction, the high spin 
fraction q at the transition temperature Tc is different from the usual value q = 1/2 
obtained for J — 0. 

In the present situation, i.e. case (3), both m and q undergo phase transitions. 
For the magnetization, when the temperature increases, magnetic property changes 

as the dia > para — > ferro-magnetic and then paramagnetic state. If J is large, 

the low temperature phase becomes ferromagnetic and then we have the change; 

the ferro > to para-magnetic states. Concerning the high spin q, its temperature 

dependence can be either a first-order transition or a continuous conversion. These 
results obtained by MFT will be studied also by MC simulations in the next section. 

To investigate the phase diagram of the SC system with magnetic interactions; we 
first expand the equations (18) and (19) around the magnetic transition temperature 
Tc (second-order case) or Tm for which m ^ 0. After some analysis, we arrive easily 
to the following equation, giving the transition line in the {zJ, T) plane 



with d = jj^, tc = ^ and j = j^- Here the reduced critical temperatures corre- 
sponding to Tq/zK and Tyi/ zK are denoted by the same quantity te- 
la order to observe the re-entrant phase transitions with the present model, the 
function j = f{tc, d) must have a minimum. Due to the self-consistent structure of 
equation (21), it is difficult to predict analytically the coordinates (^crit; Jcrit) of this 
minimum. Let us consider separately the left and right sides of the equation (21), 
and denote these two quantities by /i(tc) — and f2{tc) = tc In (2g^^^ . At the 

critical temperature, we must have fi{tc) — /2(^c) and /{(^c) — f^i^cji where the /' 
denotes the derivative of / {tc) ■ These two relations lead to the following transition 



d = tela 

J 




(21) 
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line for the critical points 

d + icrit - 1 = icrit In (2^^) • (22) 

Solving numerically the last equation, using d — D/zK — 150/90 = 5/3 and g — 10, 
we obtain jcrit ~ 0.64, which gives zJcrit — 57.6 K. This value is in very good 
agreement with the thermal behavior of the order parameters of Fig. 4. These 
results are summarized in the phase diagram of a SC system including magnetic 
interaction between the HS units, depicted in Fig. 5. A reentrant phase transition 
appears for zJcth — 57.6K < zJ < zJq — 210K where jo = 2.33. The magnetic 
phase appears from Tq to Tm- A magnetic second-order phase transition occurs at 
Tm- On the other hand, a change from a second-order to first-order magnetic phase 
transition at Tq takes place while J becomes large from Jcrit to Jq. When we choose 
larger K, T^q becomes smaller and the transition at Tq can become of first order 
(see next section). 

FIG. 5. The phase diagram of a SC system with magnetic interactions in the 
plane. The coordinate (jcritj *crit) denotes the end point of the reentrant phase transition. 
A reentrant phase transition appears for z Jcrit = 57. 6K < zJ < zJq = 210K and magnetic 
order appears from Tq to Tm. A magnetic second-order phase transition occurs at Tm- 
On the other hand, magnetic transition at Tm changes from second order to first order 
with J. 

Let us now discuss the case of small J for which m is always zero. In this case, 
although the magnetic interaction may cause a metastable magnetic phase at low 
temperatures, the system simply undergoes a transition from the diamagnetic state 
(LS) to the paramagnetic state (HS). This situation corresponds to the simple spin- 
crossover phenomena. The transition can be smooth or discontinuous as the first 
order transition. The latter is the case of the experimental results observed in PBA, 
like Nao.44Co[Fe(CN)6]o.73 • 2.7 H2O [56]. 

It is worth noting that many authors have investigated the phase diagram of 
the original BEG model [57] . Most of works on BEG model concerned the case of 
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antiferromagnetic biquadratic interaction {K < 0) [58,57,59] because it is expected 
a rich phase diagram due to competing interactions. It was found that BEG model 
exhibits various comphcate phase transitions, such as successive phase transition, 
re-entrant and double re-entrant double transition against the temperature. Here, in 
the present study, both magnetic (J) and biquadratic (K) interactions are positive 
(ferro), and as long as these interactions are "ferromagnetic", the original BEG 
model does not show re-entrant phases. However, in the present study we obtain a 
re-entrant behavior with J > and K > 0. On the other hand, the degeneracy of 
the HS state (s = ±1) combined to the ligand field D act as an effective temperature 
dependent " anisotropy" {D — k-eT In g) which is strongly reduced when temperature 
increases. As a consequence, in a rough approach, the original phase diagram (in the 
plane — j at constant D) of the original BEG model (in which g = 1) is changing 
depending on the temperature. This effect constitutes the basic mechanism of a new 
competition leading to obtain re-entrant phase transitions. 



IV. MONTE CARLO SIMULATION 

In the previous section, we have studied the structure of the phase diagram as a 
function of the lattice coupling K and magnetic couphng J in MET. In this section, 
we check the results obtained by MET using a Monte Carlo method where the effect 
of fluctuation is taken into account. Here, we adopt an extended scheme of the 
Glauber dynamics because we study only static properties of the model. We deal 
with three-dimensional system {z = 6) with 20x20x20 sites through this section. 
The dynamical properties such as the relaxation process after rapid changing of 
temperature or photo-irradiation will be reported in the separate paper [38] where 
the types of dynamics will be carefully discussed. 

Taking into account the degeneracy of the states, u and r for the low and high 
spin states, respectively, we adopt the following transition probability between the 
initial state S and the final state S': 
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where P{S) is the degeneracy of the electronic state, represented by the fictitious 
spin 5" (5* = ±1 for the HS state and 5* = for the LS state). Thus, we have 
P(0) = and P(l) = P(-l) = r. Here, E{S) is given by 



This energy denotes the energy of the corresponding site, which can be HS or LS and 
it includes elastic (trough the "biquadratic" interaction) and magnetic contributions 
due to the interactions with the neighbors. The transition probability W{S —>■ S') 
satisfies the detailed balance for the canonical distribution for the degenerate state, 
we define the following quantity m' to study magnetic order in the MC method 



where (■ ■ ■) means MC average and N is the system size (20'^). 

FIG. 6. (a) Temperature dependence of q (O) and m' (A) obtained by the Monte Carlo 
simulation, (b) Temperature and field dependence of q for the first order phase transition. 
Circles (Q) denotes the case h = 0, (A) h = 9 and (□) h = 18. (c) Temperature and field 
dependence of m' for the first order phase transition. Note that the magnetization shows 
non-monotonic behavior for nonzero h. 

In Figs. 6, we show the temperature dependence of q (Q) and m' (A). Here the 
temperature is raised up from T = 5 to 80K sequentially, and is reduced to T = 5K 
again. In Fig.6(a), smooth change of q and no change of m' is depicted for J — 0, 
h ^ 0, g ^ 10, K ^ 15K and D = 150K. The transition temperature (Teq = 35. IK) 
is larger than the critical temperature (Ti = 16. 9K) of the corresponding Ising model 
and no phase transition occurs. We find that the plotted circles for both upward 
and downward processes are well superimposed at each temperature. In Fig. 6 (b) 
and (c), we plot q and m' in the case of the first order phase transition, respectively. 



E{S) ^-JS^Sj-hS - KS'^ ^ -Sj + (L> - k^T In g) S'^. 



(24) 



j j 



m — 




(25) 
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where c, = 10, X = 28K, and D = 150K, giving Teq = 22.0K (< Ti = 31.5 K), for 
J = 0, /i = (O) and for ^ ( A, □ ) . 

We performed relatively short steps (200MCS) at each temperature for Figs 6, 
which is enough to get equilibrium for case (a), leading well visible hysteresis loops 
of q for case (b). 

The shape of the hysteresis changes with the duration time in principle. Indeed, if 
we change the temperature infinitesimally slow, the system would change in a quasi- 
static way and the hysteresis would disappear. However, in practice, the nature of 
metastability is robust and the shape of the hysteresis changes little in the case with 
5000 MCS at each temperature. In the Fig.6(b) we also confirm the field-dependence 
of the shape of the hysteresis obtained in the previous section by MET. Triangles 
(A) denote the case h = 9 and squares (□) h = 18. Shift of hysteresis curve to 
the low temperature side is observed and the shift is more prominent in downward 
process. 

FIG. 7. Temperature dependence of g(0) and m'(A) for K = 15K, D = 150K, g = 10, 
and J = 12K obtained by the Monte Carlo simulation. Inset shows the heat capacity (•). 
Note that the magnetic interaction drives the first-order SC transition. 

Here, we study the effect of magnetic interaction on the spin-crossover phase 
transition. We consider the case where the magnetic interaction J is nonzero. If we 
include weak magnetic interaction, the qualitative nature of the order parameter q 
does not change, although the interaction causes a metastable branch of magnetically 
ordered state at low temperatures. This effect will be reported when we study 
the dynamical properties of the system [38]. Here, we investigate the change of 
equilibrium phases which were studied in the previous section in MET. In order to 
change the structure of phase diagram we need a strong magnetic interaction which 
is comparable to K. In Eigs.7, temperature dependence of q and m' is depicted for 
K = 15(K) and D = 150(K), where Teq > ^i, with J = 12(K). At each temperature, 
2000 MCS are performed to obtain the data. We find that the magnetization appears 
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at intermediate temperatures. Decreasing the temperature, magnetization appears 
continuously at the critical temperature Tm — 43K where the phase transition is of 
the second order as was predicted in MFT, and vanishes discontinuously in the first 
order phase transition at Tc — 32K, where the HS fraction changes discontinuously. 

FIG. 8. Phase diagrams in the coordinate {J/K,T/zK) obtained by the Monte Carlo 
simulation. Closed squares denote Tm, closed triangles Tq, and closed circles the SC 
transitions, (a) for K = 15K, D = 150K, and g = 10 (b) for K = 22K, D = 150K, and 
5 = 10. 

Now, we study dependence of the critical temperatures Tm and Tq as a function 
of J for the case K = 15(K), D = 150(K) and g = 10, where Teq > Tj. In Fig.8(a), 
the dependence is depicted, where the closed circles denote the temperature at which 
q — 1/2. The transitions depicted by the closed circles occur smoothly. The closed 
triangles denote Tq and the close boxes Tm- The magnetic transition at Tc changes 
from second order to first order with J between J/K — 0.73 and J/K — 0.8. 

As another example, we also give the phase diagram in Fig. 8(b) for the case 
K = 22(K) and D = 150(K), where Teq < Tj. Then the closed circles denote the 
phase boundary of the first order phase transition. Technically we obtain Teq and Tc 
of the first order phase transition as follows. We prepare an initial condition where 
half of the system is in LS state (S = 0) and the other half is in the magnetic HS 
state {S — 1). We perform MC simulations with different random number sequences 
for given J and observe the relaxation of q and m'. When all samples go to LS (HS 
para, HS magnetic) phase, we regard the parameter set as belonging to the LS (HS 
para, HS magnetic) phase. We could determine the border within the precision of 
the size of the closed circles and closed triangles. These phase diagrams are very 
similar to those obtained by MFT in the previous section. It is interesting to study 
the intermediate parameter region between (a) and (b), where the line connecting 
the closed circles with the smooth transition (case (a)) changes to that with the first 
order phase transition (case (b)). 
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V. SUMMARY AND DISCUSSION 



The phase diagram of the system with the spin-crossover phase transition and 
the magnetic phase transition has been studied using an extended BEG model with 
degenerated levels by means of the mean-field theory and Monte Carlo simulations. 
This model can be a model for both SC and PBA materials in the framework of a 
common description. The present model can be also seen as an extension to 3 states 
of the Wajnflasz and Pick model for spin-transition [24]. This unified approach 
contains both "elastic" and magnetic interactions between the molecular species. 
These two kind of interactions obey to different symmetries, and then involve two 
different order parameters. Indeed, the terms responsible for the SC phenomena 
are the temperature-dependent ligand field {D — k^T In g) and the biquadratic 
—KSfSj interaction, which modulate the energy between HS and LS states, leading 
to the gradual or first-order transition. On the other hand, the magnetic interaction 
—JSiSj, assumed here as ferromagnetic, acts only between the HS units, which 
breaks the symmetry, and leads to the second-order or first-order magnetic transition 
when it competes with the SC transition. 

We found that when the magnetic interaction is weak, it does not cause a mag- 
netic phase transition at even very low temperature at equilibrium. It only causes 
a metastable magnetic state at very low temperature. This is consistent with the 
observation in PBA. The properties of the metastable state will be reported in the 
study of dynamics [38] . It has also been found that the magnetic interaction causes 
a magnetic order in the HS phase at intermediate temperatures, i.e. Tq < T < Tu- 
The magnetic order disappears at a lower temperature than Tq- 

Even in the condition of K and D that q{T) changes smoothly when J = 0, i.e., 
Teq > Ti, the magnetic interaction drives the spin crossover transition which causes 
the first-order transition for large J. 

As we have referred in the text, the dynamical properties of the present model 
are also very interesting and important in the light of recent experimental extension 
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of photo magnetization. [60-63] They would be related to a description of the photo- 
magnetic behavior of PBA, in particular, of the relaxation properties of the photo- 
induced magnetic state. There, the relaxation of two coupled magnetic and elastic 
order parameters, which have different time scales, is involved. [38] 
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